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ABSTRACT 

Superbursts are very energetic Type I X-ray bursts discovered in recent years by long term monitoring of X-ray 
bursters, and believed to be due to unstable ignition of carbon in the deep ocean of the neutron star. In this 
Letter, we follow the thermal evolution of the surface layers as they cool following the burst. The resulting 
light curves agree very well with observations for layer masses in the range 10 25 -10 25 g expected from ignition 
calculations, and for an energy release > 10 17 erg per gram during the flash. We show that at late times the 
cooling flux from the layer decays as a power law F oc f" 4 ' 3 , giving timescales for quenching of normal Type 
I bursting of weeks, in good agreement with observational limits. We show that simultaneous modelling of 
superburst lightcurves and quenching times promises to constrain both the thickness of the fuel layer and the 
energy deposited. 

Subject headings: accretion, accretion disks-X-rays:bursts-stars:neutron 



1. INTRODUCTION 

Type I X-ray bursts from accreting neutron stars in low 
mass X-ray binaries involve unstable thermonuclear burn- 
ing of accreted hydrogen (H) and helium (He) (Lewin, van 
Paradijs, & Taam 1995). In the last few years, long term mon- 
itoring of X-ray bursters by BeppoSAX and the Rossi X-Ray 
Timing Explorer (RXTE) has revealed a new class of very en- 
ergetic Type I X-ray bursts, now known as "superbursts" (see 
Strohmayer & Bildsten 2003; Kuulkers 2003 for reviews). 
The 10 42 erg energies and several hour durations of super- 
bursts are 100-1000 times greater than usual Type I bursts. In 
addition, they are rare: so far 8 have been seen from 7 sources, 
with recurrence times not well-constrained, but estimated as 
~ 1 year (Kuulkers 2002; in 't Zand et al. 2003; Wijnands 

2001) , instead of hours to days for usual Type I bursts. 

The current picture is that superbursts are due to unstable 
ignition of carbon at densities p ~ 10 8 -10 9 g cm" 3 . Hydro- 
gen and helium burn at p ~ 10 5 — 10 6 g cm" 3 via the rp-process 
(Wallace & Woosley 198 1), producing chiefly heavy elements 
beyond the iron group (including nuclei as massive as A = 104; 
Schatz et al. 2001), but with some residual carbon (mass frac- 
tion X c ~ 0.01-0.1) (Schatz et al. 2003). Cumming & Bild- 
sten (2001) (hereafter CB01) showed that this small amount 
of carbon can ignite unstably once the mass of the ash layer 
reaches ~ 10 25 g (see also Strohmayer & Brown 2002). This 
fits well with observed superburst energies for Xc ~ 0.1 and 
an energy release from the nuclear burning of 1 MeV per nu- 
cleon. The low thermal conductivity of the rp-process ashes 
gives a large temperature gradient and ignition at the required 
mass (CB01). The heavy nuclei may also photodisintegrate 
to iron group during the flash, enhancing the nuclear energy 
release (Schatz, Bildsten, & Cumming 2003). Therefore su- 
perbursts offer an opportunity to study the rp-process ashes. 

Previous authors used one-zone models to estimate the 
time-dependence of the flash (CB01; Strohmayer & Brown 

2002) . In this paper, we present the first multi-zone models 
of the cooling phase of superbursts. Unlike normal Type I 
bursts, the time to burn the fuel is much less than the convec- 
tive turnover time. We therefore assume that the fuel burns 
locally and instantaneously in place, without significant verti- 



cal mixing. We do not calculate ignition conditions, but rather 
treat the amount of energy deposited and the thickness of the 
fuel layer as free parameters 2 . In §2, we describe our calcu- 
lations of the subsequent thermal evolution of the layer, and 
present a simple analytic model which helps to understand the 
numerical results. At late times, the flux evolves as a power 
law in time rather than the exponential decay found by CB01 
for a one-zone model. In §3, we use the long term flux evolu- 
tion of the layer to predict the timescale of quenching of Type 
I bursts after the superburst, and compare to observations. 

2. TIME EVOLUTION OF THE SUPERBURST 

After the fuel burns, the cooling of the layer is described by 
the entropy equation 
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where the heat flux is F = -K(dT / dr), and t v is the neutrino 
energy loss rate. The layer remains in hydrostatic balance, 
in which case a useful independent coordinate is the column 
depth y into the star (units: g cm" 2 ), where dy = -pdr, giving 
a pressure P = gy. The surface gravity is g = (GM/R 2 )(l+z), 
where 1 +z = (1 -2GM /Rc 1 ) -1 ! 1 is the gravitational redshift 
factor. In this paper, we assume M = 1.4M and R = 10km, 
giving z = 0.31 and gu = g/10 14 cm s" 2 = 2.45. 

To find the temperature profile just after the fuel burns, we 
deposit an energy E mc = E\i 10 17 erg g" 1 throughout the layer. 
Since carbon burning to iron gives w 10 18 erg g" 1 , we expect 
En = 1 to correspond to Xc ~ 0.05-0.1, depending on how 
much energy is contributed by photodisintegration (Schatz 
et al. 2003). At each depth, we calculate the temperature 
of the layer, Tj, from J^ 1 cpdT = E DUC , where 7] the initial 
temperature. Following CB01, a simple analytic estimate of 
Tf is as follows. The electrons are degenerate and relativis- 

tic for p > 10 7 g cm" 3 , giving pY e = 1.5 x 10 8 g cm" 3 P^ 4 
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Fermi energy Ep = 2.7 MeV P 2 g , and pressure scale height 

2 This is a similar approach to Eichler & Cheng (1989) who studied the 
thermal response of a neutron star to energy deposition at different depths. 
However, the transient events they consider are less energetic than super- 
bursts. 
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FIG. 1 . — Temperature and flux profiles for yt = 10 12 g cm 4 , £17 = 1, and 
* = 10 minutes (solid lines), 1 hour (dotted lines), 1 day (short-dashed lines), 
and 10 days (long-dashed lines) after ignition. In the lower panel, the t = 1 
hour (1 day) flux profile is shown scaled by a factor of 30 (300). 
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FIG. 2. — Flux from the surface against time for different values for the 
energy release £"17 and layer thickness yi, (in g cm -2 ). 
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FIG. 3. — Cumulative energy released from the surface as a function of 
time, for yt = 10 12 (solid lines) and 10 13 gcirT 2 (dotted lines) and in each 
case (bottom to top) £17 = 1, 2, and 3. 



H = y/p = 61m P l 2 ^Y e /g XA , where P = P 26 1 26 erg cm" 3 . The 
heat capacity cp is determined mainly by the electrons, cp w 

ir 2 (Y e k B /m p )(k B T/E F ) = 2.6 x 10 7 erg g" 1 K" 1 T 9 Y e P^ /4 . In- 
tegrating and assuming 7} <C Tf, we find the convectively sta- 
ble temperature profile Tf = 3.6 x 10 9 K E\!j 2 P^Y e 1|/8 , insen- 
sitive to depth and depending mainly on £17. 

Our thermal evolution code uses the method of lines, in 
which the right hand side of equation (1) is differenced over 
a spatial grid, and the resulting set of ordinary differen- 
tial equations integrated using a stiff integrator. We choose 
a uniform grid in sm\\~ x {\ogy /yb), which concentrates grid 
points around the base of the layer y/,, resolving the initial 
temperature discontinuity. We place the outer boundary at 
y = 10 8 g cm" 2 , and set flux oc T 4 there; at the inner boundary, 
typically y w 10 14 g cm" 2 , we assume vanishing flux. We as- 
sume the layer is heated before the flash by a 10 21 erg cm" 2 s" 1 
flux from the crust, and that all the fuel burns to 56 Fe. Our 
results are not sensitive to the details of the grid, or bound- 
ary conditions (for times longer than the thermal time at the 
top zone). We calculate the equation of state, opacity, neu- 
trino emissivity, and heat capacity as described by Schatz et 
al. (2003). 

Figure 1 shows temperature and flux profiles 10 minutes, 1 
hour, 1 day, and 10 days after ignition for a model with En = 1 
and yb = 10 12 g cm" 2 . Figure 2 shows a series of lightcurves 
for different yb and £17. At early times, as the outer parts of 
the layer thermally adjust, the radiative flux depends mostly 
on £17. At late times, after the cooling wave reaches the base 
of the layer, the flux depends mostly on yb, and falls off as 
a power law F oc f" 4 / 3 . Figure 3 shows the cumulative en- 
ergy release fory^, = 10 12 and 10 13 g cm" 2 and £17 = 1-3. In 
the first few hours, the energy released from the surface is 
w 10 42 ergs, the exact value being mainly sensitive to £17, 
rather than depth. A significant fraction of the heat is ini- 
tially conducted inwards and released on a longer timescale, 
as pointed out by Strohmayer & Brown (2002). 

The physical reason for the late-time power law flux decay 
is that as time evolves, the peak of the temperature profile 
moves to greater depths where the thermal timescale to the 
surface is longer 3 (see Fig. 1). The simplest analytic model 
is a slab with constant thermal diffusivity D, whose tempera- 
ture is perturbed close to the surface, for example at a depth 
x = a (where x = is the surface). For a delta-function per- 
turbation initially, the temperature evolution is given by the 
Green's function 4 

, x 1 • , / ax \ ( x 2 + a 2 \ 

T(x, t) = ■ sinh - — exp (2) 

V^Dt \2Dt) \ 4Dt J 

and the surface flux is F oc (dT/dx) x=0 oc r 3 / 2 exp(-T/f), 
where r is the thermal time at the initial heating depth 
r = 4a 2 /D. For an initial "top hat" temperature pro- 
file T(x < a) = 1, T(x > a) = 0, the surface flux is F oc 
{r/t) 1 ! 2 [l -exp(-r/f)] . For t < t, before the cooling wave 
reaches the base of the layer, F oc f" 1 / 2 ; for t > t, the solution 

3 A similar problem is ohmic decay of crustal magnetic fields, where power 
law decay is also expected (Sang & Chanmugam 1987; Urpin, Chanmugam, 
& Sang 1994) 

4 A simple way to obtain this result is to apply the method of images to 
the Green's function for an unbounded domain T(x,t) oc r 1 / 2 exp (-x 2 /4Dt) . 
Eichler & Cheng (1989) derive a similar result for a power law dependence 
of conductivity on depth, which also shows self-similar behavior at late times 
(see Lyubarsky, Eichler, & Thompson 2002 for a recent application to cooling 
of SGR 1900+14 after an outburst). 
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is independent of the initial temperature profile, and F oc t ~ 3 ^ 2 . 

The numerical results show a similar behavior, although 
with different power law indices. The relevant timescale in 
this case is f c00 ] = H 2 /D, where D = K / pep. The electron con- 
ductivity is K = n 2 n e I^T fom+iic, where = Ef/c 2 , and v c 
is the electron collision frequency. When electrons dominate 
the heat capacity, the thermal diffusivity takes the particu- 
larly simple temperature-independent form D = c 2 /3v c . For 

electron-ion collisions, v c = 9.3 x 10 16 s _1 '{Z 2 /A)A ei /Y e 
(e.g. see Appendix of Schatz et al. 1999), giving 
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(see also eq. [10] of CB01), where we insert the appropriate 
numbers for 56 Fe composition. The simple "top hat" solution 
for constant conductivity motivates a fit to the numerical so- 
lutions, 



F 25 = 0.2 t, 
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where tht = t/l hour. For t > f 000 i, £25 = 0.13 E^ 2 (t /f coo i) 4/3 - 
The transition from F oc f" 02 to F a ?" 4 / 3 occurs when 

f/ ^cool ~ ^17 

Equation @ fits the numerical results to better than a factor 
of two for models without substantial neutrino emission. As 
emphasised by Strohmayer & Brown (2002), neutrino cool- 
ing is important for large carbon fractions: it depresses the 
flux at t « 5-10 hours for the models with En = 2 and 3, 
y = 10 13 g cm" 2 in Figure 2. Whenever neutrinos dominate 
the cooling, the peak temperature is large enough that emis- 
sion is by pair annhilation. A good fit to the neutrino energy 

loss rate is e u w 10 4 erg g" 1 s~ l Tg 2 y l2 ^ 2 , giving a cooling time 
?„ = cpT/e u = 2.5 x 10 12 s y 5 J, 4 T 9 ~ 10 . Inserting the peak tem- 
perature from equation (2) gives t v sa 300 hrs E^. Neutrinos 
dominate when t v < t coo \, or when £17 > 2.3 y^ 20 - 

3. COMPARISON TO OBSERVATIONS 

The cooling curves in Figure |2] compare well with ob- 
served lightcurves, including a rapid initial decay on hour 
timescales, followed by an extended tail of emission (as ob- 
served following some superbursts, e.g. KS 1731-260, Ku- 
ulkers et al. 2002; Ser X-l, Cornelisse et al. 2002). We wiU 
present a detailed comparison with the observed superburst 
lightcurves in a future paper. The initial decay from the peak 
depends mostly on £17, and so it should be possible to con- 
strain the amount of fuel consumed in the superburst. Our 
models do not resolve the peak itself, since this depends on 
the details of how the burning propagates out to the surface; 
however, for En > 2-3, the flux exceeds the Eddington flux, 
^Edd = 3 x 10 25 erg cm" 2 s _1 /(l +X), where X is the H frac- 
tion, for timescales of minutes. Superburst peak luminosities 
are generally less than the Eddington luminosity (Kuulkers 
2003), implying £17 < 2. The one exception is the super- 
burst from 4U 1820-30, which showed dramatic photospheric 
radius expansion lasting for several minutes (Strohmayer & 
Brown 2002). This is consistent with the proposal that this 
source, which accretes and burns He rich material, produces 
large quantities of carbon (Strohmayer & Brown 2002; Cum- 
ming 2003a). The transition to the late-time power law oc- 
curs after t « 4 h E^ Al y\ 2 A (eq. [5]), which corresponds to 
£25 ~ 0.13£ 2 7. It may therefore be possible to measure the 
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FIG. 4. — Observed and predicted quenching timescales of normal Type 
I bursts following a superburst. We show the predicted quenching time for 
yi, = 10' 2 (solid lines) and 10 13 g cirr 2 (dotted lines) as a function of both 
the critical flux needed to stabilize H/He burning F cl j t , and the accretion rate 
(given in terms of F cl ; t by eq. |5j). For each value of yi,, we show curves 
for (bottom to top) £17 = 1, 2 and 3. The observed upper or lower limits on 
tq Uenc h, and the estimated M, are taken from Table 1 of Kuulkers 2003, with 
an updated value for4U 1636-53 (Kuulkers, private communication). 



power law decay using superburst tails, although this depends 
upon being able to subtract out the underlying accretion lumi- 
nosity, Faccr,25 ~ 0.1 (M/OTMndd), in a reliable way. 

Another way to probe the late-time cooling is to use the 
remarkable observation that Type I bursts disappear (are 
"quenched") for f que nch ~ weeks following the superburst 
(e.g. Kuulkers 2003). CB01 proposed that the cooling flux 
from the superburst temporarily stabilizes the H/He burning. 
An estimate of the critical stabilizing flux, Fdu is as follows. 
The condition for temperature fluctuations to grow, and un- 
stable He ignition to occur, is vtT, a = r/e coo i (Fushiki & Lamb 
1987), where e^a is the triple alpha (3a) energy production 
rate, e coo i is a local approximation to the cooling rate, and v 
and 77 are the respective temperature sensitivities. For a large 
flux from below, the He burns stably before reaching this ig- 
nition condition, at a depth where the time to accumulate the 
layer equals the He burning time, y/m = YQ^aje^a, where Y 
is the He mass fraction, m is the local accretion rate per unit 
area, and Qj, a = 5.84 x 10 17 erg g" 1 = 0.606 MeV per nucleon 
is the 3a energy release. At the transition from unstable to 
stable burning, both criteria are satisfied at the base of the 
H/He layer. Using the first condition to eliminate e?, a from 
the second, and writing e coo i w F/y, gives F = lymQ^Y /r] = 
6.2 x 10 22 erg cm -2 s- 1 (»V'«Edd)(>'/0.3)(^/4r7). Some of this 
flux is provided by hot CNO burning of accreted H, Fh ~ 
e H y = 5.8 x 10 21 erg cm" 2 s^ysCZ/O.Ol) (Cumming & Bild- 
sten 2000; Z is the metallicity); the remainder is £ cr ; t = F-Fh- 
This estimate agrees well with a more detailed calculation us- 
ing the ignition models of Cumming & Bildsten (2000), in 
which we find F alX ps mQ^ a w 0.7 MeV per accreted nucleon, 
almost independent of M. Therefore, 

£crit,22 ~ 6 (»i/»iEdd) (5) 

(see also Paczynski 1983a; Bildsten 1995). Equation (5) in 
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the limit t ^> t eoo \ gives 

ao, 77-3/4^3/8 , , 3/4 _-3/4 ,-,3/8 , £ s 

fquench = JO 'cool / 'crit.22 i 17 = " da y S ^12 ^crit,22 fc 17 ' W 

which gives f quen ch in terms of the thickness of the layer and 
energy release. 

Figure 5 compares the predicted and observed quenching 
times. The observations of f que nch an d accretion rates (used to 
find Feat from eq. 0) are taken from Kuulkers (2003) (except 
for 4U 1636-53, which has a revised upper limit of 23 days, 
Kuulkers private communication). The observations are up- 
per or lower limits only: nonetheless, the general agreement 
is very good and supports the quenching picture suggested 
by CB01. There is much to learn from a careful compari- 
son of superburst lightcurves and the corresponding quench- 
ing times, separately constraining both En and yj. 

4. SUMMARY AND CONCLUSIONS 

We have presented the first multi-zone models of the cool- 
ing phase of superbursts. The flux decay is not exponential, 
but power-law (eq. [5]). For t < f coo i, where f COO ] is the cool- 
ing time at the base of the layer, the flux depends mostly 
on the energy release En, and is insensitive to depth: the 
inwards travelling cooling wave does not yet know that the 
layer has a finite thickness. For t > f coo i, the flux decays as 
a power law F oc f" 4 / 3 , independent of the initial temperature 
profile. The power law decay at late times gives predicted 
Type I burst quenching times of weeks (eq. [6]), consistent 
with observational limits. Future comparisons of both super- 
burst lightcurves and quenching times with observations will 
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constrain both the thickness of the fuel layer and the energy 
deposited, particularly when combined with models of normal 
Type I bursts from the same source (Cumming 2003a,b). 

There is still much to be done in terms of theory. Perhaps 
the most important issues are the physics of the rise (which 
sets the initial condition for our simulations), and production 
of the fuel. Important clues to the first are the observed pre- 
cursors to superbursts, which may be normal Type I bursts ig- 
nited by the superburst. Recent progress has been made on the 
second, with indications from both theory (Schatz et al. 2003; 
Woosley et al. 2003) and observations (in 't Zand et al. 2003) 
that stable burning may be required to produce enough car- 
bon to power superbursts. A self-consistent model of H/He 
burning, followed by accumulation and ignition of the ashes 
may require a better understanding of the transition from un- 
stable to stable burning observed in normal Type I bursting 
(e.g. Comelisse et al. 2003). 
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